Chern numbers hiding in time-of-flight images 
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We present a technique for detecting topological invariants - Chern numbers - from time-of-flight images of 
ultra-cold atoms. We show that the Chern numbers of integer quantum Hall states of lattice fermions leave their 
fingerprints in the atoms' momentum distribution. We analytically demonstrate that the number of local maxima 
in the momentum distribution is equal to the Chern number in two limiting cases, for large hopping anisotropy 
and in the continuum limit. In addition, our numerical simulations beyond these two limits show that these local 
maxima persist for a range of parameters. Thus, an everyday observable in cold atom experiments can serve as 
a useful tool to characterize and visualize quantum states with non-trivial topology. 



I. INTRODUCTION 

Topological insulators (TI) and superfluids are many-body 
quantum systems with energy gaps in the bulk but topolog- 
ically protected, gapless excitations on the boundary ||Tl|2l- 
Each class of topological phases is characterized by a topolog- 
ical invariant which can be either an integer or a binary quan- 
tity. For example, the integer quantum Hall insulators 1 3 1 are 
distinguished from other two-dimensional insulators with the 
same symmetry by the Chern number, an integer that capture 
the global properties of all occupied bands |4|. The Chern 
number coincides with the number of chiral edge modes, and 
the Hall conductance in units of the conductance quantum. A 
central problem in exploring new forms of topological matter 
is how to extract the topological invariants from experiments, 
and unambiguously determine the topological phase. 

Experiments with ultracold atoms has brought us opportu- 
nities to creating topological phases of matter in parameter 
regimes unreachable in solids. In addition to band-structure 
engineering using optical lattices and tuning the atom-atom 
interaction by Feshbach resonance, artificial magnetic fields 
for neutral atoms have been produced by either rotation Q or 
atom-laser coupling f6^ 'T\, making way to realizing quantum 
Hall states of cold atoms |8|. Moreover, non-Abelian syn- 
thetic gauge fields |7 9| and spin-orbit coupling flOl have 
been demonstrated. These developments motivated theories 
on properties of ultracold topological matter, e.g., the frac- 
tional quantum Hall effect of interacting bosons 151 [TTl ITSl . 
the quantum spin Hall effect of fermions 19J, and the quantum 
anomalous Hall states of p-orbital fermions lfT3l . 

While cold atom systems offer new detection techniques 
and great tunability, it is challenging to measure the topologi- 
cal invariants directly. Measurements of local density of states 
(of the edge states) or mass transport are marred with com- 
plexity in contrast to solid state systems. Thus, creative meth- 
ods have to be invented. For example, Ref. |8| proposed to 
extract the quantized Hall conductance from the derivatives 
of the in-situ atom density distribution by the Stfeda formula. 
Here, we show that the Chern number of a lattice-quantum 
Hall system of fermions can be simply determined from the 
structures (tipples) in the momentum distribution n(k) as 
measured in time-of-flight (TOE) images for a range of pa- 
rameters. Our results transform n(k), a simple everyday ob- 



servable of cold-atom experiments, into an unexpected, use- 
ful tool for characterizing and visualizing topological states 
of matter 

Figure [T^ summarizes the main results of this paper ob- 
tained for the Hofstadter model (14] . We highlight the pa- 
rameter regimes where the structure in n(k) is most evident: 
in the continuum limit (Fig.[T]3) or for highly anisotropic lat- 
tices (Fig.[T];). In particular, we provide a transparent under- 
standing of n(k) for large hopping anisotropy. In this limit, 
states with energy at the band edges dominate, and they are 
dimerized (their localized wave functions are peaked at pairs 
of spatially separated lattice sites). The size of these "Chern- 
dimers," in units of the lattice spacing d, is equal to the Chern 
number of the corresponding gap. As a result, n(k) varies si- 
nusoidally (Fig. [T];) with a period determined by the Chern 
number. Furthermore, our detailed numerical studies ver- 
ify that such distinctive features in the momentum distribu- 
tion persist beyond the Chern-dimer and the continuum limit 
(Fig.[2j; and[2|l). However, the coiTespondence between the 
number of local maxima in n(k) and the Chern number is not 
always precise for arbitrary parameters. 



II. ULTRACOLD FERMIONS IN SYNTHETIC MAGNETIC 

FIELD 

We consider the Hofstadter model llT4l of spinless 
fermionic atoms on a square optical lattice in the e^-ey plane 

H = -Y^ is,s' {e''-'flfs' + h.c.) -Y^fisflfs- (1) 

Here, is the creation operator for a fermion at site 
s — {sx,Sy}; the nearest-neighbor hopping strengths = 
ts.s+e^, and ty ~ ts^s+sy are real and positive; ds^s' — 

e A ■ dl depends on the vector potential A 11211 for parti- 
cles with charge e; and /ig is a local chemical potential. With 
cold atoms, the hopping anisotropy A — ty/tx can be tuned 
by varying the intensity of the lasers giving rise to the optical 
lattice potential |6|. For an applied magnetic field B = Bbz 
normal to the 2D plane, a = d^eB/h is the magnetic flux 
per plaquette in units of the flux quantum. We study homo- 
geneous gases with /ig /i, and also finite systems in the 
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presence of an isotropic harmonic trap potential, captured by 
Ms = Mo ~ muP'<Ps^ /2. 

Artificial magnetic fields have been implemented for 
bosons 1 15 1, and it is straightforward to apply these techniques 
to fermions. With their smaller fine-structure splittings, the 
alkali fermions (^Li and ^°K) are more susceptible to heat- 
ing from spontaneous emission than some bosons (^^Rb and 
^■^■^Cs), but for the h/tx^y ^ 240 /is timescales ||22ll associated 
with the major energy gaps in the spectrum (see Fig. [2^), the 
expected w 200 ms lifetime in '"^K should be acceptable; see 
Ref (9 1 for an alternate solution. 

We first summarize the well known properties of the ho- 
mogeneous Hofstadter model [Eq. ([T]l] with commensurate 
flux a = p/q, where the integers p, q are relatively prime. In 
the Landau gauge |23| suitable for current experiments |6|, 
the vector potential A = axGy scales the unit cell along 
by a factor of q to a magnetic unit cell; this also creates a 
magnetic Brillouin zone (MBZ) with G (— 7r/g, Tr/q] and 
ky E (— TT, 7r]. The Hofstadter model is diagonal in the single- 
particle basis 'i's-,,Sy = e^^''^''''~^^''^''''^s^{ky), provided ^pg^ 
satisfies El 

gjfc.^r _^^_^_^-^K^r _^+2Xcos{2nSxa+ky)ipl = -Er^Pl 

(2) 

with energy in units of t^- The index r — 1,2, ...g labels 
linearly independent solutions, and V^s^+g = V'l^- The single 
energy band at a splits into q sub-bands, and the energy 
spectrum has q — I gaps 1 14]. When the chemical potential /i 
is inside the j'-th energy gap, the system is an integer quantum 
Hall state with Chern number c^, defined by the Thouless- 
Kohmoto-Nightingale-Nijs (TKNN) formula 101 and related 
to the quantized Hall resistance. 

Usually we consider TOF images as approximate represen- 
tations of the momentum distribution n(k) obtained by imag- 
ing the 2D atomic column density after a t^oF period of free 
expansion. In most cases, the final spatial coordinate of an 
atom is X = v^tof where v is the velocity |16|, however, 
vector potentials complicate the interpretation of TOF images 
because the velocity is v = /ik/m — eA/m. In Ref. ifTTll . we 
showed that suddenly turning the gauge field to zero is equiv- 
alent to a brief electric field that takes the initial momentum 
fik and maps it to a new final velocity v' = ^k/m. As a 
result, detected TOF images do approximate momentum dis- 
tributions [24.1 . The crystal momentum distribution 

n(k)=^(/t/,,)e'''-(-^'), (3) 

s,s' 

can be reconstructed from TOF images using the recipe in 
Ref [18|. Here, the average (• • •) is over the ground state. We 
also define a ID crystal momentum distribution by integrating 
71 (k) along ey, 

n{kx) = ( dkyn(k). (4) 

J —7T 

For homogeneous systems, we computed n{kx) by numer- 
ically diagonalizing Eq. (j2]) using a sufficiently large unit cell, 
while for finite trapped systems we diagonalized Eq. ([T]l. Our 



numerical results are complemented by analytic analysis for 
low flux and large anisotropy. In what follows, we discuss 
in turn the behaviors of n{kx) in different parameter regimes, 
starting with an isotropic lattice and low flux. 

III. OSCILLATIONS IN MOMENTUM DISTRIBUTION 

Figures [TJ) and|2j3 show n{kx) in this isotropic low flux 
limit, at a = 0.05 and A 1. As the chemical potential is in- 
creased up to fi — 0, successively higher Hofstadter bands are 
filled, with Chern numbers Cr — r where r — 1, 2, 3... labels 
the highest filled band. Correspondingly, n{kx) has exactly c,. 
local peaks on top of a smooth envelope, i.e., counting peaks 
in TOF images unambiguously determines the Chern number, 
and therefore which quantum Hall state the system is in. To 
understand this, we recall that for a ^ 1, the cyclotron orbit 
is large compared to the lattice spacing, so the system resem- 
bles a free Fermi gas in magnetic field, and the narrow bands 
in the Hofstadter spectrum map to Landau levels. For r filled 
Landau levels, n{kx) is the summed square modulus of the 
momentum-space wave functions, 

r-l 

n(fc,) cx ^(2''^.!)-liy2(fc^£5)e-('=^^«)^ (5) 

where H^, are the Hermite polynomials, and £b = \fhfeB is 
the magnetic length. Due to the node structure of Hy, n{kx) 
has exactly r local maxima located at the zeros of iJr+i- This 
is reflected in our numerical results as plotted in Fig. [T|5 for 
the lattice system. As is depicted by the shaded regions in 
Fig. [T^, this association is valid in the entire low flux regime 
and is insensitive to A. 

At larger a the Landau level picture breaks down, but 
the correspondence between the number of local peaks and 
the Chern number appears to persist for a range of parame- 
ters. Fig.|2j; and |2jl show examples of n{kj:) for large flux 
a = 0.45 and A = 1 where well resolved peaks do remain. 
More generally, as a increases (for /i residing within the same 
energy gap), the distinctive features of n{kj.) are generally re- 
tained, but compared to the low flux regime, the local peaks 
(except at k^ — 0) move to higher momenta and can be sub- 
sumed into the background structure. In some instances, these 
peaks turn into shoulder-like structures. 

Intriguingly, as the hopping anisotropy becomes large, 
small wiggles on a background turn into regular and pro- 
nounced oscillations regardless of a. Fig. |2j; and |2]l illustrate 
the progression of n{kx) from A = 1 to A = 10, showing 
pronounced oscillations develop with increasing A. Most im- 
portantly, the simple rule reappears in the limit of A 1: the 
number of oscillations in n{kx), i.e., the number of peaks for 
fcj. e [— TT, tt), is again exactly equal to the Chern number c^. 
As is emphasized in Fig. [I];, n{kx) asymptotically approaches 

n{kx) cx r ± ^ cos{crkx) (6) 

for chemical potential inside the r-th gap. The one-to-one cor- 
respondence, or mapping, between the number of oscillations 
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and the Chem number for anisotropic lattices (Fig.[T^: Chern- 
dimer regime) is one of our key observations which we now 
explain analytically. 



IV. CHERN-DIMERS 

The A ^ 1 limit is tractable in perturbation theory lfT9l 
where the solutions to Eq. (|2]) are simple. When = (i.e., 
A — ?> oo), there are q bands each with dispersion Er{ky) = 
—2ty cos{27r ar+ky ), r G {1, . . . , g} (see Fig. [5] where p — 1 
and q = 5). These eigenstates are extended along Gy and lo- 
calized to single sites as a function of x. As defined in Eq. Q, 
the r-th band describes states tbl = Sg r localized at site 
Sx = r and displaced in ky, as identified in Fig.p] The bands 
intersect at numerous degeneracy points k* for states r — 
and r' — s'^. The spatial separation along e^, of these lo- 
calized states is £ = s'^ — s^- These band-crossings occur 
at q — 1 discrete energies (right axes of Fig. [3]l with exactly q 
crossings having the same £ at each energy. For finite A, hy- 
bridizes the states at each crossing, leading toq{q — l) avoided 
crossings, and {q — 1) gaps in the band structure. Near each 
degeneracy point the system is well described by the effective 
Hamiltonian L19il . 

i?2(k) = Ae[cos{kx£)a-x+sm{kxi)cry]+ve{ky-ky)az. (7) 

d'x.y,z are the Pauli matrices in the basis of localized 
states \r) and |r'); and the A; -independent constants are 

A. = (-i)^n'^z.i^,{i [Erik;) + Er>ik;)]~En{k;)}-' 

and Vi = 2A [sin(27rar + fcp - sin(27rar' + k*)]. The 
eigenstates are "Chern-dimers", superpositions of |r) and \r') 
with wave function 

^f^ = ^{S,^^, + e'P^ds^y), (8) 

where /3+ = —£{kx + n) + tt and /3_ = —£{kx + tt) are rel- 
ative phases for the upper and lower band edges, respectively. 
Starting from Eq. (j7]l it is straightforward to show |19| that 
the Chern number of each band gap is = £. Thus, for suffi- 
ciently anisotropic hopping, dimerized states form at the band 
edges with spatial extent along e^; equal to the Chem number 
of the corresponding gap. 

The formation of dimers of size ^ = at the edges of the 
r-th gap implies a "hidden spatial correlation" described by 
correlation function C{1) = X]s(/s /s+/e^), which is peaked 
at I = £. For chemical potential in the r-th gap, C{1) asymp- 
totes to a delta function C{1) — >■ S{1 — £), as A — ^ oo. Be- 
cause /3_|_ — /3_ = TT for each crossing, the net contribution to 
C{1) from dimerized states associated with gaps fully below 
the Fermi energy is zero. Since C{1) and n{kx) are a Fourier 
transform pair, n{kx) = y\ C(Z) cos(fca;Z)/-\/iV; this yields 
the cos{crkx) term in Eq. (j6]). The r offset in Eq. (j6|l comes 
from the C{1 = Q) contribution. 



V. EFFECT OF TEMPERATURE AND TRAP 

Our numerical simulation based on Eq. (1) confirms that 
the structures in n(fcj.) persist, and remain visible, for trapped 
systems at finite temperatures. Due to the finite band gap of 
the quantum Hall insulator, thermal fluctuations only slightly 
smear the peak (or oscillation) structure in n{kx) as long as 
T is small compared to the gap of interest, which is of order 
t/2 for the major gaps in the Hofstadter spectrum. Fig. |4] 
shows examples of n{kx) for realistic trap and temperature 
parameters. Thus, our proposal of detecting the Chern number 
is experimentally feasible. 



VI. Q-MATRIX AND MOMENTUM DISTRIBUTION 

Finally, we point out a general connection between n(k) 
and the Q-matrix ll20l . which plays a fundamental role in 
the classification of topological insulators and in turn is re- 
lated to the Chern number of a quantum Hall insulator For a 
band insulator with band structure -E, (k) (with energy mea- 
sured from the chemical potential), the spectral projection op- 
erators -P+(k) and P-(k) project onto states above or below 
the chemical potential, respectively. The Q-operator Q(k) = 
P_i_(k) — P_(k) describes a "flattened" Hamiltonian that is 
smoothly deformed from the original ll20l . The momentum 
distribution is related to the spectral projection operator of the 
filled bands p(k) = _P_(k). In matrix form, p is related to Q 
by 

Qab(k)^Sab-2pabi'k). (9) 

As an example, for the Hofstadter model at flux a — p/q, 
Pab{k) = (/i(k)/fc(k)), where a,b e {1, . . . , g} are the site 
indices within the magnetic unit cell, and k is the crystal mo- 
mentum within the MBZ. n(k) as defined in Eq. ^ is sim- 
ply n(k) = Ea6Pa6(k). By Eq j9|, for given ky = k*y, 
the behavior of h{kx,k*) as function of k^ is determined 
by the Q-matrix along cut ky ~ k* in the BZ. So is the 
ID momentum distribution n{kx). On the other hand, the 
Chern number is nothing but the winding number of the off- 
diagonal block of the Q-matrix in chiral basis, also known 
as the transition function, along a cut that splits the BZ |20|. 
Thus, on general grounds, phase winding responsible for 
nonzero Chern number implies nontrivial structures in the 
momentum distribution function. Explicitly, for the asymp- 
totically exact effective Hamiltonian /f2(k), the Q-matrix 
is Q(k) = {vjkl + A2)-i/2ij2(k). The Chern number 
Cr — £, for the off-diagonal block of Q-matrix is simply 
^ 6*'^=^. Meanwhile, from the Q-matrix, we obtain n(kx) — 
q [1 — sign(Af ) cos(kx£y\- This proves analytically our previ- 
ous result Eq. (j6]l. The oscillations adiabatically evolve into 
local peaks for the isotropic lattice, as explicitly shown in 
Fig.|2}; and|2}J. 

In contrast with Chern numbers, n(kx) is not a topologi- 
cally invariant quantity. Still, the clear signatures of Chern 
numbers in TOE images that we find here makes momentum 
distributions into an unexpected tool for exploring topologi- 
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cal states of matter. In addition to edge states and Hall con- 
ductance, ripples in the momentum distribution and dimerized 
states of anisotropic lattices whose spatial extent encodes the 
topological quantum numbers provide a new way to charac- 
terize and visualize non-trivial topological states. Our study 
opens the possibility of using momentum distributions to ex- 



plore other topological states of matter including topological 
superfluids. 
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(a) Schematic map of parameter space 




(b) Isotropic, iiearly-coiitiiiuumm: X = 1, a = 0.05 
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FIG. 1: (Color) Relation between the Chern number Cr and structure 
in TOF data as a function of magnetic flux per unit cell a in units 
of the flux quantum h/e, and hopping anisotropy A. (a) Schematic 
map of the a-A plane showing the Chern-dimer and the continuum 
regimes. Numerical results for the remainder of parameter space will 
be shown in Fig. [2]and[3] (b,c) Normalized ID crystal momentum 
distributions n{kx), see Eq. (|4]l, computed at three different chemical 
potentials corresponding to Cr — 2, 3, 4; in each case, the number of 
peaks or oscillations is equal to Cr- The data in (b) is for an isotropic 
A = 1 lattice at a = 0.05, and the data in (c) is for an anisotropic 
A = 1000 lattice at a = 0.45. ka: is measured in units of 1/d, where 
d is the lattice spacing. 
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FIG. 2: (Color) ID crystal momentum distributions n{kx). (a) En- 
ergy spectrum of Eqs. ([T]) and ^ for A = 2. (b) Evolution of n(fca;) 
with increasing chemical potential fi al X = 1 and a = 0.05 show- 
ing Cr increasing in tandem with the number of oscillations in n{kx). 
(c,d) Evolution of n{kx) with increasing anisotropy A at a = 0.45 
for (c) Cr = 2 and (d) Cr = 3. These curves depict the distributions 
approaching the cosine function with increasing A. The unit for kx 
is 1/d. 




FIG. 3: (Color) Energies in the anisotropic limit for p = 1, 
g = 5 showing the perfect sinusoidal energies when A — >■ oo 
(solid black curves), and the avoided crossings when A = 10 (solid 
blue curves). The hybridized states at these crossings correspond 
to "Chern-dimers," and the red-circled example has spatial extent 
5sx = 1. The unit for ky is 1/d. 
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FIG. 4: (Color) The momentum distribution function n[kx) for the 
isotropic lattice at temperature T — 0.2t and a trap potential with 
muP'iff' = 10~'^t. /io at the trap center is set to — 0.75t, — 0.5t for 
a — 0.2, 0.3 respectively. The twin peaks, indicated by the arrows, 
allow to determine the Chern number near the trap center to be 2. 
The unit for fc^^ is 1/ d. 



